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Abstract 

In the preceding paper (arXiv:1103.0329 [math-ph]) we treated the Jaynes-Cummings 
model with dissipation and gave an approximate solution to the master equation for the 
density operator under the general setting by making use of the Zassenhaus expansion. 

However, to obtain a compact form of the approximate solution (which is in general 

complicated infinite series) is very hard when an initial condition is given. To overcome 
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this difficulty we develop another approach and obtain a compact approximate solution 
when some initial condition is given. 

This paper is a sequel to the preceding one pQ. In the paper we treat the Jaynes- 
Cummings model with dissipation (or the quantum damped Jaynes-Cummings model in 
our terminology) once more and study the structure of general solution from a mathematical 
point of view. 

We want to apply it to Quantum Computation and Quantum Control which are our final 
target [2], [3]. As a general introduction to these topics see for example [3] and [5]. 

We expect that our study will become a starting point to study more sophisticated models 
with dissipation in a near future. 

Let us start with the following phenomenological master equation for the density operator 
of the atom-cavity system in jS] : 

J^p = —i[H, p] + ji ^apa) — -(a)ap + pa) a) j- + u |a^pa — -(aa^p + paa)) j (1) 

where H (for simplicity we write H not Hjc in [TJ) is the well-known Jaynes-Cummings 
Hamiltonian (see [7]) given by 



— cr 3 <g> 1 + u l2 <g) a)i 
^ + u N fia 



H = -^0"3 <8> 1 + u l 2 <S> a) a + fi (cr+ <g> a + er_ <g> a 1 ") 



fiat _^i +Wo AT 



(2) 



with 




\ / 1 \ / 1 

= 1 I • = I ' a3 = ' 12 = 

10/ \ -1 / \ 1 

, and a and at are the annihilation and creation operators of an electro-magnetic field mode 
in a cavity, iV = a) a is the number operator, and p and v (p > v > 0) are some constants 
depending on it (for example, a damping rate of the cavity mode). 

Note that the density operator p is in M(2; C) ® M(J") = M(2; Mp 7 )), namely 

p = I P0 ° P01 I £M(2;I(J)) (3) 
P10 P11 



where M(J-) is the set of all operators on the Fock space J 7 defined by 



JF = Vect c {|0),|l),|2),|3),...} 

oo oo 

^Cn\n) | ^|c n | 2 < oo J>; |n) = ^i=r|0) 



. n=0 n=0 

and 1 in 02]) is the identity operator. 

Now, we decompose (j2J into diagonal and off-diagonal parts 

H = H d + H off =\ z + ! fi) 

-tf + tu N J \ fiat 

and rewrite (jTJ as 

— p = —i[H d , p]+p |apa^ — -{a^ap + pa^a) |+z/ |a^pa — -{aa)p + paa^) | — i[H ff, p\. (5) 
Namely, the main part is 

—i[H d ,p] + p ^apa^ — -((Jap + pa^a)! + v pa — -(aa^p + paa)) 
, while a kind of perturbed one is 

-i[H off ,p\. 

This is the main difference between [T] and this paper. Although the form may be not 

standard it is very useful to calculate when some initial conditions are given. 
Let us write down the equation (JS]) in a component-wise manner. Then 

Poo = -iu) {Np 00 - p 00 N) + p japooa 1 - ^(a f ap o + Poo^a) j + v |a f pooa - ^(aatpoo + Pooaa f ) 
-ifl(apio - poia 1 ), 

Poi = -iua(poi + Npoi - poiN) + p |apoia t - ^(a 1 apoi + PoiaM j + v ja 1 p 01 a - ^(aaVoi + Poiaa t ) 
-iQ(apu - pooa), 

Pio = -iua{-pio + Np w - p w N) + p ^ap Q1 a^ - ^{a)ap m + p 01 a^a)^ + v ja 1 p 01 a - i(aaVoi + Poiaa^) 
-ifl(a f poo - pno)), 

pn = -iu} (Np n - pxxN) + p japna 1 - ^(a^ap n + pne^a) j + v ja f pii<2 - ]^(ao) p lx + pnaa r ) 

-^(aVoi - Pioa) (6) 



where py = (d/dt)pij as usual. 
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Here we use some technique used in [T] (see also [8] and [9] ) , which is very useful in some 
case. For a matrix X = (x^) G M(JF) 

( \ 

^00 ^01 ^02 ' ' ' 



X 



X\Q X\\ Xyi • • • 

X20 #21 ^22 " " " 

V : ; '''J 

we correspond to the vector X e J r dim c .F ag 



X — [Xij) — > X — (xqo, 3^01, ^02, • • • ; xio, xix, x±2, ■ ■ • ; X20, x%i, X22, 



(7) 



where T means the transpose. Then the following formula 



EXF = {E® F T )X 



holds for E,F,X e M{JF). 
This and equations (jH]) give 



Poo 



1 



-iu (iV®l-l®iV)+//<{a® {a r f - -(a ] a ® 1 + 1 ® a ] a) \ + 



?<■ <; a ] ® a T - ^(aa f <g> 1 + 1 <g> aa ] ] 



poo - in (a ® lpio - 1 ® (a f ) T poi) 



P01 



1 



-2W (l®l + iV®:L-l®AO+//-ja® (a 1 ") 1 " - -{a ] a ® 1 + 1 ® a) a) \ + 



v \ a} ® a T - -(aa f ® 1 + 1 ® aa f ) 



P01 — ifi (a ® lpu - 1 ® a T poo) 



P10 



1 



-iu (-1 <g>l + iV<g}l-l<g}iV)+p<Ja<g> (a ] ) T - - (a f a ® 1 + 1 ® a ] a) \ + 



v \ a ] <g> a T - -(aa ] ® 1 + 1 ® aa f ) 



P10 - ifi (a 1 " ® lpoo - 1 ® (a f ) T Pii) 



P11 



1 



-iu (N ® 1 - 1 ® N) + p |a ® (a 1 ) 7, - -(a f a ® 1 + 1 ® a f a) }> -f 
z/ I a f ® a T - ^(aa f ® 1 + 1 ® aa ] ) 



pn - iVl (a f ® lpoi - 1 ® a T pio) (9) 
because 1 and N = a^a are diagonal (1 T = 1, N T = N). 
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Next, in order to rewrite matrix elements in terms of Lie algebraic notations used in [8] 
we set 



a f <ga T , K_ = a®(a ] ) T , K 3 = ~(N <g> 1 + 1 ® N + 1 <g> 1), 



K + 

K = N <g) 1 - 1 <g) N. 



(10) 



Then it is easy to see 



[K 3 ,K + ]=K + , [K 3 ,K_] = -K_, [K + ,K_] = -2K 3 , 
[K ,K + ] = [K ,K-] = [K ,K 3 ] = 0. 



(11) 



Namely, {K + , K_, K 3 } are generators of the Lie algebra s«(l,l), see for example [TUj as a 
general introduction. 
If we set from (131) 



Poo Poi 



Pio Pu 



Poi 
Pio 



(12) 



we obtain the following "canonical" form 
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at 



p = (x + r) P 



(13) 



with 



X 



-iu K + L 




/ 



y = -in 





-1 ®a T 
a f <g> 1 






-iojQ — iuj K + L 

iuo — iuoKo + L 



-1® (a)) T o®l 



a f <g> 1 -1 (g) a T 







a <g> 1 
-1 ® (at) T 








-iu K + L 



\ 



and 



L /i \ a <g> (a ] ) T - -(a'a <8> 1 + 1 <8> a ] a) 1 + v \a ] <g> a T - ^(aa 1 g) 1 + 1 <g> aa f ) 



fjL — V 



+ uK + + - (/x + i/)AT 3 



where the fundamental relation aa^ = a^a + 1 = iV + 1 and a simplified notation w in place 
of WqI <g) 1 have been used. 

Let us note once more that X is not anti-hermitian, while K is anti-hermitian. 
Since the general solution of the equation (|T3|) is given by 



p(t) = e'( x+y )/)(0) 



(14) 



in a formal way we must calculate the term e^ x+Y \ which is in general not easy (to obtain a 
compact form is almost impossible). For that the following Zassenhaus formula is convenient. 



Zassenhaus Formula We have an expansion 

t(A+B) _ 



e -^{2[[A,B],B] + [[A,B],A\} e ^[A,B] e tB e tA ^ 



(15) 



The formula is a bit different from that of [TT 
In this paper we use the approximation 



e t(X+Y) ^ e ^[X,Y] e tY e tX_ 



(16) 



Let us calculate each term explicitly. 
[I] First, we calculate e tx . The result is 

/ e t(-iLu K +L) Q 

e -iu t e t(-iu) K +L) 


\ 



JX 








e iu}ot e t(-ioJoKo+L) 









\ 



Q e t(-iu) K +L) 



J 

and fortunately in [8] the term ^(-^oKo+l) ^ as Deen calculated exactly. Namely, 



(17) 



,t(-^oK +L) _ „^t„G(t)K + -iu tK -2log(F{t))K 3p E(t)K- 



where 

-^sinh(^t) 

u—v \ 2 ) 



E(t) 



cosh(^t)+^sinh(^) 



F(t) = coshl^tU^sinhf^ 



2 J /i — v 
^sinh(^t) 

= cosh(^)+^^h(^)- (19) 
This is a kind of disentangling formula, see for example [TO] as a general introduction. 
If from (HHD 

= gt(-»o;o^o+i)^o) = e i£ i il * e G(t)^+ e -^o^o-21og(F(t))X3 e ii;(t)X-^Q^ ^0) 

then the original form is given by 



r W = WE^r( at ) ri { ex p({-^- 1 °g^w)}iv)x 

^ ' n=0 

^ ^Z! a - r (0)(at) m i exp ({icu t - log(F(t))}N)}a n . (21) 



, m=0 



See [§]. Next, we list some results from [S] for the latter convenience, 
(i) If r(0) = |0)(0| then 



r(t) = e jQe^ N . (22) 



(ii) If t(0) = where |a) (a G C) is a coherent state defined by 

| a j = g^-^lo) = e -^e aat |0) a|a) = a|a)) 

then 

r(t) = (l-G(t))eH 2e - ( ^ )tlo ^W e I J. (23) 

The main part (which corresponds to the classical one) 

ae -(^^) W + - e -(^-^ )t a 

appears in the formula, see Appendix. This derivation is not easy, so see |5] for further 
details. 
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[II] Second, we must calculate e tY . We decompose Y into two parts 

Y = -in{Y x - Y 2 ) 

where 



/ 



a (8)1 

o®l 

a f ®l 

a)® I 



\ 



/ 



Y? 



1 <g> (at) T 
1 ® a T 




V 



1 ® a T 

From (J2I|), fl25l) and [1] it is easy to see 

a 






1 ® (a f ) T 




\ 



(24) 



(25) 













a 1 









exp — zfit 



at 



exp I zfit 



a 



at 



cos(fitv aat 



t= sin (fit "\/ aat) a 



i—rr= sinffitVataW cos (QtV at a) 

VaTa v ' v ' 

cos(fitv^aat) sin (fit V aat) a 

v V aaT v ' 

t= sin (fit \/ ata)at cos (fit V at a) 

'at a v ' v ' 



T 



(26) 



t 

[III] Third, we must calculate e~^ x,Y \ For the purpose we first calculate [X, Y], which is 
relatively easy. Note that Yi and Y 2 commute from (|25|) . 
Then the result is 

[X,Y] = -in{[X,Y 1 ]-[X,Y 2 ] 
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where 



and 



[X, YV 



' A ^ 



A 
B 
5 



( C ^ 



D 
c 
D 



C 



= -va} <S> 1 + 



z/1 (g) a T , 



5 



M; ^ t 
— z — a 1 



(a f ) T , D = na®l 



1 d 



1 <g> a T . 



It is easy to see that 



[A,C] = [A,D] = 0, [B,C] = [B,D] = 0, 



so we can conclude that [X, Yi] and [X, Y 2 ] commute. 
Since 

e §[xy\ = e -4m,Yi] e 4m,Y2] 

we can calculate each term easily. The result is 

„-*#fi[X,Yi] _ 



/ 





sin(4^V A B]4)-B 





cos(^fiVAB) 


sin(4^\/SA)B 



/BA 



and 



/AB 







/AB 



cos(^-Q\/iL4) 









cos(^-fi-/BA) / 



( cos(^-flVCD) 






cos(^-nVDC) 










/CD 








IDC 
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In last, we shall restore the result to original form. For the purpose, ignoring the term 

e Vl x ' y ] we set 

Pit) = e tY e tx p(0) = e tY h (*) , h (0 = e tx p(0) (27) 

and 

' Poo(0) Poi(0) 
Pio(0) Pn(0) 



P(0) 

Then from [I] pi{t) becomes 



where 



~Pi(t) = | (12) | (28) 

(21) (22) 



^ ' n=0 

£ ^!" ampoo(0)(atr I 6XP ({ ' Wot " lo sW))}^)K, 

,m=0 m ' J 

(12) = e ^*^^_y_( a tHexp({-< Wb t-tog(F(0)}N)x 

^ ' n=0 

^a m p O i(0)(at) m exp ({iwtf - log(F(t))}iV)K, 



,m=0 



(21) = e-^^^(at)«{exp({-^ t-log(F(t))}iV)x 

^ ' n=0 



V ^pa m Pl0 (0)(at) m 1 exp ({iu t - \og(F(t))}N)}a n , 



to! 

,m=0 
^-" 1 00 



( 22 ) = ^TYE^r( at n ex P({-^-log(F(t))}iV)x 

^ ' n=0 

£ _Lj_ a - pil(0 )( a t)- exp ({iu, t - log(F(i))}A0K 



to! 

. m=0 



and from [II] p{t) becomes 



cos(QWN + 1) -?^=r sin(fitV-W + l)a 
-i^ sin(fitViV)a t cos(fitViV) 

cos(fW-W + 1) ^T^f sm(ttWN + l)a 
sin(fit v / iV)a t cos(My/N) 
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or by making a slight modification in terms of af(N) = f(N + l)a 

cos(£lty/N + 1) -i^w+i sm{QWN + 1 



a 



pit) =| ; ■ ■ - v*n * • ■ < ]h{t)x 

cos (QWN + 1) ia^-sm(nty/N) 



i a 



= sm{nty/N + 1) cos{QtVN) 



(29) 



By making use of this formula let us calculate an important example. The initial state is 



Example 

1 l |0)(0| 
P(0) = - [ 

\a){a 



1 / |0> <0| 

P(0 =o ( 3 °) 



where |a) is a coherent state in [I]. 
Then the result is 



~Pi(t) = U A I (3D 



2 



B 



where 



U— V . 



-t 

A = e 2 C \QR G(t)N 



F(t) 

(see (122D and (El) and 



G(t)<ae U J a'+ae V 2^ u Ja-N 



2 

where 



«.) = U (11) <12) I (32) 



(21) (22) 



11) = cosjtttVN + 1) A cos^tViV + 1) + = ^ sin(OtyiV + l)oBa t = = sin(QtyiV + 1), 

12) = icos(nty/N + l)Aa-^=sm(nty/N) - i - == = sm(Qty/N + l)aB cos(QtVN), 



(21) = -i-j— sin (QtVN)a! A cosjntV N + 1 ) + i cosjntVN) R^t 1 = sin (fltV /V + 1 ) , 
y N v N + 1 



(22) = — L sin(fitv / iV)a t Aa-^ sin(fitv / iV) + cos(QtVN)B cos(fitv / iV). 
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These forms are compact and comparatively beautiful. Though we can of course calculate 
another example we stop here. 



In this paper we reconsidered the Jaynes-Cummings model with dissipation from a dif- 
ferent point of view and constructed a compact approximate solution when some initial 
condition was given. It is very fresh as far as we know. We will leave a further construction 
to readers who are interested in this topic. As for the preceding works see [12], [13] and |14j . 

mi- 

We conclude this paper by stating some future prospects. Our real target is the following 
master equation : 

— p = —i{Hn, p] + p |apa^ — -{a)ap + pa) a) \ JrV | a V a ~~ ~^ (aa) p + paa)) 

where Hr is the Rabi Hamiltonian (without RWA (Rotating Wave Approximation)) given 
by 

H R = <8> 1 + W0I2 ® + Qai ® (a + at) 

f + u N Q(a + at) 
ft(a + at) -f + u N 

We call this the Rabi model with dissipation. The Jaynes-Cummings model (which 
is an approximate model with RWA) has some weak points (see for example [16J and its 
references), so we must treat a more realistic model like this. In the following paper (s) we 
will attack this model. 



Appendix 

In this appendix we review the solution of classical damped harmonic oscillator, which 
is important to understand the text. See any textbook on Mathematical Physics. 
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The differential equation is given by 

x + 72 + u?x = (7 > 0) (33) 

where x = x(t), x = dx/dt and the mass is set to 1 for simplicity. In the following we treat 
only the case u > 7/2 (the case u = 7/2 may be interesting). 
Solutions (with complex form) are well-known to be 

x± (t) = e -(^V^W)\ 
so the general solution is given by 

x(t) = {ae-(^^^y + ae^V^^y} X (0) 

= {« e -(i + ^v /IZ ^)* + ae-(^^v / M^)*}x(0) (34) 

where a is any complex number. 

If , ~f/2u is small enough we have an approximate solution 

x(t) « {ae-^ +iu) > + ae-^-^H x(0). (35) 
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